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The usage of vortex labyrinths fields and Talbot lattices as optical dipole traps for neutral atoms is 
considered for the wavelength of trapping radiation in the range 0.98 -r 2.79 /im. The square vortex 
lattices generated in high Fresnel number solid-state microchip lasers are studied as a possible 
realization. The distribution of light field is obtained via dynamical model based on Maxwell- 
Bloch equations for class-B laser, discrete Fox-Lee map with relaxation of inversion and static 
model based on superposition of copropagating Gaussian beams. The spatial patterns observed 
experimentally and obtained numerically are interpreted as nonlinear superposition of vortices with 
^3 helicoidal phase dislocations. Such patterns are approximated analytically by a sum of array of 

vortex lines and additional parabolic subtrap. The separable optical trapping potential is proposed. 
The factorization of macroscopic wavefunction have led to analytical solution of Gross-Pitaevski 
5— l , equation for ensemble of quantum particles trapped in vortex labyrinth formed by spatially - periodic 

array of Laguerre-Gaussian beams. 
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PACS numbers: 03.75.Lm., 42.50.Vk., 42.60.Jf, 42.65.Sf.,47.32.Cc 



I. INTRODUCTION 



Optical dipole traps for neutral atoms [l[ are a subject of considerable interest. A certain amount of proposals of 
increasing complexity of trapping EM-field appeared recently. Apart from a relatively simple geometrical patterns like 
i— — 1 1 standing plane- wave gratings 0, evanescent wave mirrors @, toroidal traps which utilize the intensity distribution of 
Laguerre-Gaussian (LG) beams in beam waist 0] a several proposals were made which utilize the arrays of Gaussian 
beams, both phase-locked [1, @] and unlocked ones 0- The interference inherent to phase- locking provides multiply 
connected configurations of intensity distribution, phase gradients and electromagnetic (EM) momentum density (allj. 
The winding EM momentum distribution is the cause of the angular momentum transfer to macroscopic bodies (e.g. 
■ dielectric ball) [l(| and to trapped BEC as well [ill ]. 

The EM field configuration under consideration is based upon properties of self- imaging optical fields [l2l [i~3Tj . The 
difference between phase-locked array of zero-order Gaussian beams and phase-locked array of optical vortices (OV) 
obtained experimentally in the near field of solid-state microchip laser [l4| is that the latter consists of array of 
parallel vortex lines with opposite circulations and topological charges Iem (TC) (fig. [T]) []|. In contrast to earlier 
proposals where individual loading and addressing of trapping sites was considered this OV array have sophisticated 
^Zh . configuration of intensity [5J and EM momentum density of trapping field. Due to this configuration the BEC trapped 
in EM field might have macroscopic wavefunction of complex form composed of array of superfluid vortices (SFV) . 
The cause of SFV formation is light-induced torque experienced by isolated resonant atom which interacts with 
Laguerre-Gaussian beam having phase singularity. It was show by Allen et al.Q the value of torque in saturation 
limit is T = H£em T- The origin of torque is due to nonzero azymuthal EM momentum component. The azimuthal 
Doppler shift corresponding to such motion had been onserved[9j- This torque might have an appreciable value even 
in nonresonant case, although it is significantly reduced by the multiple of (A/r) 2 , where A - is detuning, T - is 
linewidth[8|. As a result an azimuthal component of EM momentum is transferred to atom in such a way that it 
will move around phase singularity - the direction of rotation is fixed by TC of trapping beam. This is the cause of 
circular motion of BEC confined in isolated toroidal trap. 

In current paper the maximally transparent and simple analytical description presented which describes the transfer 
of angular momentum from optical vortex array (OVA) to trapped BEC. As a first step in section 2 the distribution of 
light intensity and phase in the near field of solid-state microchip laser are computed. The configuration of intensity 
and phase of EM-field of OVA has complicated, multinode structure. In sections 3 and section 4 the spatially periodic 
arrays of Gaussian and Laguerre-Gaussian beams of the first order are compared. Next in section 5 in order to get 
analytic solution for macroscopic wavefunction 'F trapped by optical vortex array, the special optical pancake-like 
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potential (fig. [TJ will be constructed. In this trapping potential the separation of variables in Gross-Pitaevsky 
equation^ becomes possible. Next the procedure applied previously to elongated sech 2 - profile optical trap [15| will 
be used. Due to special ajustment of potential it is possible to capture an analytic approximation for \1/ in the form of 
superposition of elementary equispaced vortices. The topological charges and angular momenta of adjaiccnt vortices 
are counter-directed in contrast to "rotating bucket" trap where angular momenta of SFV in bucket are co-directed 
[ltjl . H3] • The numerical modeling via split-step FFT technique will be performed to check this approximate formula. 
The obtained spatial distribution of W modulus and argument will be used for analysis of classical field of velocities of 
a trapped atom. The complex form of constructed macroscopic wavefunction W trapped by optical vortex array might 
be interesting from the point view of diminishing the decoherence induced by environment in topological quantum 
computing 



II. SQUARE VORTEX LATTICES 



Recent advances in controlling the dynamics of solid-state microchip lasers [3] offer the possibility of reliable control 
of spatiotemporal optical patterns. Compared to semiconductor lasers the host medium composed of dielectric crystal 
doped by neodimium Nd +3 or other rare earth ions (Er +3 , Tm +3 , Ho +3 , Yb +3 ) have smaller gain and smaller self- 
phase modulation. Because of smaller nonresonant losses the heating of host crystal by radiation is not dramatic. The 
changes of the geometry and birefrigence of host crystals, curvature and reflectivity of output couplers and spatial 
distribution of optical pumping give efficient control over mode structure. The square vortex lattices (SVL) observed 
in quasi plane parallel cavity are of special interest. These lattices demonstrate high degree of spatial coherence 
: the relaxation oscillations of class-B high Fresnel number solid-state microchip laser [3| with Fresnel number in 
the range Nf w 100 -j- 1000 are characterized by single peak at frequency about (VTir c ) _1 ( 0.3 -j- 1.2 MHz ).This 
is the firm evidence of single-longitudinal and single transverse mode behaviour. The former theoretical analysis was 
based on parabolic equation, resulted from adiabatic elimination of polarization from standard set of Maxwell-Bloch 
equations [19| : 

dE{r,t) E{r,t) ic = acN Q L a E(f,t)(l + i8uT 2 ) 

dt t c 2k ± [r ' ' 2L r (l + aT lC e \E\ 2 /tkj) {) 

where r c - is photon lifetime in cavity, k=2ir/X - wavenumber, a - stimulated emission crosssection, Suj - detuning 
of lasing frequency from center of gain line, iVo - density of inverted resonant atoms per unit volume, T\ - inversion 
lifetime (longitudinal relaxation lifetime), L a - thickness of active medium, L r - length of resonator, c - speed of 
light, £o - dielectric constant, Aj_=V^ 2 [20]. When finite gain lincwidth T 2 ~ l is taken into account in the framework 
of Swift-Hohcnbcrg equation: 



dE(f,t) , E(r,t) ir c c 
at (T c + T 2 j 2k(T c + T 2 ) 
T2 2 ,T e c. , . ^s 2cV ^ (TcN L a E(r,t)(l + iSuT 2 ) 



T C (T C + 1 2 ) Z k 2L r (1 + al\ c eo\E\ z /nw) 

, the square vortex lattices were obtained numerically [2l|. The additional term with the square of transverse laplacian 
is responsible for transverse mode selection due to finite gain linewidth.An alternative model with discrete time step 
equal to r = 2L r /c (time of bouncing of radiation between mirrors) utilizes the standard rate equations of class-B 
laser written at n-th step for electric field (H. I22I]: 

r M tor (*\ vL a N n (r)E n (r)(l + i6uT 2 ) 

E n +i{r) = f(E n {r)) = + E n {r), (3) 

inversion: 

N n+1 (r) = N n (f) + { N o(r)- Nn(rl) _ q ^ \ E ^ /nu] ^r (4) 

Ti c 

and nonlocal integral mapping evaluating the field via fast Fourrier transform at each timestep: 

/>oo P 00 

E n+l (r)= / K(r-r')f(E n (P))d 2 r', (5) 
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The kernel K for the nearly plane-parallel Fabry-Perot cavity of microchip laser with transverse filtering via aperture 
D{?) has the form [13] : 



ikD(r') , ,_ , 
K(f- r') = - \ ' exp [ik(f- r') 2 /2L r ] 
2irE r 



(6) 



The following parameters were choosen in numerical simulation: 7\ = 2-10~ 4 sec, L r — 1mm, for A r d +3 -doped 
crystals a = (1.2 — 0.6)T0 _2O cto 2 , N = 10 16 cto~ 3 , SujT 2 = 0.1. In numerical evaluation of cq. ([5]) via split-step 
FFT method [23j the mesh size in X, Y plane was 512 • 512 points. The "guard bands ratio" [26| was set equal to 
8. The main part of the field ^ was located inside central part of a mesh of 64 • 64 size - the "image area". The 
tolerance of the energy spillover was kept within ei = 0.001. The "windowing" in wavenumber space after FFT at 
each timestep was performed by usage of " Fermi-Dirac" smoothed step function (23| . The dissipation inherent to 
"windowing" in split-step FFT method was compensated by nonlinear gain. The initial conditions for field E n were 
taken as multimode random field [23| . 

Just near lasing threshold the radiation mode has a distribution of intensity as rectangular grating of bright and 
dark spots : the latter are vortex cores. The nonlocal integral mapping proved to be succesful in computation of the 
near field distribution as well. Quite unexpectably in the most runs the parallel vortex lines were obtained (fig. [J) @ 
rather than periodic array of bright and dark spots, typical to Talbot phenomenon [l2j]. The origin of parallel vortex 
lines is interpreted as nonlinear superposition of vortices with helicoidal phase dislocations. In the next section 3 and 
section 4 the possibility of approximation of this nonlinear optical vortex lattice by linear equivalent will be outlined. 



III. TALBOT LATTICES 

Consider the phase-locked rectangular lattice of zero-order Gaussian beams located at sites fj X ,jy [3 separated 
by period p, where jx,jy - are dicrctc indices corresponding x and y coordinate of a given site. Let us assume for 
simplicity that polarization is linear, thus spin of light is zero. At z — plane the electric field E is given by expression: 

E(r,0) =EoYl ex P Hr - r Jxjy \ 2 /(2d 2 )} (7) 

After paraxial propagation of distance z the electric field E(r, z) is transformed into: 

% exp [ ikz] ^ \r-fjx,j V \ 2 1 , , 

The constructive interference between ajaicent beams produces periodic interference patterns in diffrenent z-spaced 
planes. The initial periodical pattern is reproduced at so-called Talbot-distances z t — 2p 2 /A. In the intermediate 
planes z t /4m the coarser lattices with periods p/m are produced [23l [24|. As a result a 3D lattice of bright spots is 
formed (fig. 02). 

Each spot could serve as potential well for neutral atoms [l], @, d, Q , because intensity gradient will attract or repulse 
the atomic dipoles depending the sign of detuning of radiation frequency from resonance. At low frequencies (red- 
detuning) atomic dipole oscillates in-phase with trapping field and tries to align parallel to electric field. Thus potential 
energy of dipole U = —p- E(f) /2 is lower in local maxima of intensity and atoms are collected at bright spots. On 
the other hand at frequencies above the resonance (blue-detuning) atomic dipole oscillates out-of-phase and it has 
tendency to align anti-parallel to electric field. In this blue-detuned case the potential energy of dipole is higher in 
local maxima of intensity and atoms are repelled into dark regions. 

Such geometry of trapping in integer and fractional Talbot planes based on superposition of co-propagating zeroth- 
order Gaussian beams was considered earlier [BJ including possibility of manipulation of optical lattice geometry via 
mutual polarization [|| of constituting beams. 

IV. ARTIFICIAL VORTEX LABYRINTH 

Consider now the periodic array of Laguerre-Gaussian vortex beams with helicoidal phase dislocations (fig. [3]) : 

E(r, z = 0) w E Q i\r- f jxJy \)exp (- |r - f jxJy \ 2 /d 2 ) 

exp{- \r\ 2 /D 2 ) exp[i£ E MArg(r-f rx j y ) + i ir(jx + jy)}. (9) 
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The topological charge £em is assumed to be unity, the neighbouring beams (components of the sum ([9])) are ir- 
shifted. The apodization function exp (— \r\ 2 / D 2 ) is added in order to suppress the maximum of interference pattern 
at the edge of array. The beams centers are placed in the centers of rectangular grid Tj X ,j y of period p whose axes 
are parallel to X,Y. The overlapping beams produce interference pattern formed by two arrays of bright and dark 
spots rotated at 45° angle with respect to initial array of LG beams. The dark spots are of two kinds: one lattice of 
spots coincides with lattice of initial vortices Q, the other one is produced by interference and it is shifted at distance 
p/V2 along diagonal of initial lattice. The resulting interferogram has apparent 45° tilt compared to lattice of initial 
vortices (fig. [3]). 

The topological charges of dark spots (vortices) flip from one site to another. The interesting feature of this inter- 
ference pattern is the distribution of angular momentum [l3j . The initial array of LG beams carries unit circulation 
and corresponding angular momentum at each site. Such chessboard — like interferogram (fig. [3]) contains additional 
array of vortices with alternative charges. The net angular momentum tends to be close to zero, because at the central 
part of array each positively directed TC I em is compensated by four adjaicent negative ones having charge — Iem 
. The elementary cell of such a lattice consists of two 7r-shifted initial vortices with co-directed charges Iem located 
at diagonal and two counter-directed charges —Iem (7r-shifted too) placed at the other diagonal of cell. The period 
p of initial pattern in numerical recipes [j| was taken as 30/zm while width of each initial beam was set wlg = 22/im 
in order to provide significant mutual overlapping of vortices . 

The longitudinal distribution of intensity of optical vortex array is composed of periodically spaced (with period 
p) hollow tubes - vortex cores(fig. [3]). Apart from Talbot gratings which are reproduced by diffraction at z t — 2p 2 /X 
planes with corresponding period division in between planes the vortex array under consideration keeps its shape 
within Rayleigh range, i.e. at distances z < D 2 /X. Within this distance the diameter of cores is kept constant along 
z by virtue of interference of adjaicent vortices, whose helical wavefronts are perfectly matched within elementary cell 
(fig.®. 

Consider now the interaction of individual atoms with single optical beam carrying topological charge. The gradient 
force will attract the "red" -detuned atomic dipole to the intensity maximum of an isolated first order Laguerre- 
Gaussian beam - "doughnut" beam (fig. [4J, i.e. to the ring-shaped area around phase singularity. As a result the 
cloud of cold atoms will be accumulated near the maximum of intensity or "doughnut" as follows from variational 
solution of Gross-Pitaevsky equation (GPE) Q: 

„ dV(r,i) h 2 A , N ,^ ^ N 

ih y ' = — - — A&(r,t) + V ext (r)V{r,t) + 
at 2m 

+ f(r,t)|*(r,t)| 2 , (10) 

m 

with trapping potential V ext of the form 

2 9 

V ext (f) = 1 Re[a(w)] \E \ 2 [r 2 exp (-r 2 /{2w LG 2 ))] (11) 

where a(uj) is polarizability of atom[l[, m - mass of particle, a - scattering length. 

The approximate variational solution of GPE for ground state macroscopic wavefunction is Laguerre-Gaussian 
function [ill: 



(r, t) = \r\ exp (~\r\ 2 /{2w r 2 ) - z 2 /(2w 2 ) + i l BEC 0), (12) 

where Ibec is the topological charge of this vortex \1/ solution. Hence the probability distribution |\I/(r, t) \ 2 of finding 
atom some near the point r = (x, y, z) at the moment t is similar to the intensity distribution of trapping field as 
pointed out in (llj . More information could be obtained from the study of phase structure of wavefunction and 
comparison of topological charges of trapping field and BEC vortex. It is easy to realize using Madelung transform 

*(f,t) = Vp(M, 4>,t) exp (i0(H,<M)) (13) 
that there exists flow of "probability fluid" with velocity v proportional (parallel) to the phase gradient lines: 

v(r,t) = —V6(r,t) (14) 
m 

Such "flow of wavefunction" occurs around the z-axis (beam axis). The "flow" described by (fT2")) is potential and 
conservative as it should be for superfluid. This picture is complicated by vorticity of EM-momentum, inherent to 
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LG-beams [9J . Strictly speaking the rotation of classical dipole around core is accelerated by nonconservative torque 
induced by azimuthal force circulating around vortex core. For two-level atom the value of torque T is given by 
following expression obtained by Babiker, Power and Allen [![: 



T = hl *« T h + I L +*P* ] (15) 

where II - normalized trapping filed intesity, A = to — luq - detuning of trapping field frequency u) from resonant 
frequency ujq of dipole. In saturation limit T is simply : T — H£em^- In nonresonant case this torque reduces 
as T 2 /A 2 . At trapping frequencies of Nd +3 laser A = \{) 1A Hz and GaAs laser A = 10 13 Hz the D-line dublett of 
85 Rb with linewidth T — 5-10 6 Hz [l| the torque is reduced by factor depending on laser intensity II .Thus due to 
trade-off between the saturation and detuning the torque exists both near and at large detuning from the resonance. 
Consequently the loop integral of the azymuthal force over circular trajectory around core 4 (Ft x r) ■ dl is nonzero 
(fig-®- 

In classical picture because of this torque the dipole placed in the doughnut beam will move " upwards" the helicoidal 
phase staircase, i.e. it will rotate around LG-beam axis. The direction of rotation is determined by topological charge 
Iem of the trapping beam (fig. rj|. Qualitatively classical dipole is pushed by azimutally periodic electric field - it 
happens because the phase of electric field oscillations at each point in circle around the center of beam is shifted 
with respect to neighbouring points. Thus classical dipole placed at maximum of intensity of GL-beam will feel the 
" plane- wave-like" nonresonant pressure of light field. The associated azimuthal Doppler shift of the moving atom 
was observed experimentally recently [25[ . Outside the resonance the origin of torque T is interpreted in similar way: 
the Pointing vector has component S(r)t tangent to helix and local flux of electromagnetic momentum pushes dipole 
along phase gradient, i.e. in azimuthal direction. The local density of EM-momentum g = S/c 2 is proportional to 
components of Pointing vector Q: 

£qu£em c 2 2 

S(r) t = \E(r)\ 

r 

S(r) z =e c\E(r)\ 2 , (16) 

where S(r) t - is tangential component of Pointing vector, S(f) z - is axial one, |(r)| - is the distance from optical vortex 
core, lo - is frequency of trapping field. The existence of tangential component of Pointing vector became visible when 
macroscopic dielectric ball with radius larger than LG-beam core and comparable to "doughnut" radius is placed in 
center of LG-beam. The rotation of particle is induced via such torque and corresponding transfer of the angular 
momentum (Io| . 

When loaded in optical vortex lattice (fig. [5])the atomic dipole will move around the adjaicent vortex core with ac- 
celeration. The radius of rotation will increase until the dipole will approach the separatrice of the velocity fields. Next 
after certain amount of rotations around the vortex core it could jump to another vortex using bright areas between 
vortices as bridge (fig. 2]). This qualitative picture is complicated by azimuthally inhomogeneous distribution of 
intensity around each core. 

Our aim now is to show that just described classical motion have quantum mechanical counterpart. The transfer of 
angular momentum results in the specific form of macroscopic wavefunction ^ maintaining coherence all over trapping 
array: 

f(r) w (l f ~ r JX ,jy\)exp (- \f - r jxJy \ 2 /d 2 ) 
exp(- \r\ 2 jD 2 ) exp [i £ B EcArg(r - r jxJy ) + i ir(jx + jy)]. (17) 



Next section presents the method of solution of Gross-Pitaevsky equation with vortex trapping field E(r) in the 
form ©. 



V. SEPARABLE VORTEX ARRAY POTENTIAL FOR BEC 



In order to get closed form solution for macroscopic ^ it is worth to mention that azimuthal accelerating force has 
very small value, falling as A -2 under detuning from resonance Next let us introduce optical potential V ext as 
a square modulus of trapping field -E(rj_, z)[l[ . The torque T will be taken into account as a "selection" rule for 
choosing distribution of topological charges I b ec in resulting solution. 
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It was shown recently [l5| that Gross-Pitaevsky equation^ in 3D-configuration: 

d$f(r t) h 2 

. h ^] = _ n ^ + y ffyp t) + 

ot 2m 

+ ^*(f,t)|*(r,t)| 2 , (18) 

m 

admits the application of standard method of separation of variables widely used for solution of linear partial differ- 
ential equation, e.g. in quantum mechanics. The separation of variables means that wavefunction is factorized: 

*(f,i) = ^(fx,i)f||(^<) (19) 

provided the Hamiltonian is the sum of two components. First component depends on a longitudinal variable z 
and a second component depends upon transverse variables rj_ . 

Following to [la ] in order to separate variables and factorize the wavefunction let us choose trapping potential in 
the following form, as a sum of components depending on longitudinal coordinate z and transverse coordinates r± 
separately: 

2 2 

V ext (rl, z) = V z +V ± = mtJ * 2 - Re[a(cj)} \E(r ± )\ 2 + 

+ m^Y±)\ 2 (20) 

where a{ui) is atomic polarizability[l| : 

a(u) = 67re c 3 -^ ^7 3 , — jt^t (21) 

The a(u>) is assumed to be real due to large negative ("red") detuning from atomic resonance at frequency 
luq = 2.4-10 15 Hz. The imaginary part of denominator under detuning A = lu — luq for 85 Rb atoms trapped by EM-field 
at A — 1.06 4- 0.808 fim is (3-f-4.7)T0 8 times smaller than real part, so the permittivity of atom a(oj) is real with 
good accuracy [l[. The trapping field E(f]_) is assumed to be periodic function of transverse variables r± = (x,y), 
composed of LG beams placed at the nodes (jx,jy) of rectangular grid of period p © . 

The additional parabolic well with frequency uj±_ is introduced in (|20|) to get analytical solution for one trapping 
vortex. In order to avoid the interference between different trapping beams the usage of different carrier frequencies 

2 2 

is recommended for longitudinal parabolic well — "| z , vortex array beam E(r) and parabolic subtrap in ([20|) .The 
characteristic scales of potential in longitudinal direction and transverse direction are choosen to form "pancake trap" 
: u> z >> uj± (fig. [1]). The opposite case of elongated trap with ui z << u>x and "solitonic" longitudinal potential 
V z « sech 2 (z) was considered earlier using analogous procedure (la) . 

The longitudinal part of wavefunction ^u(z,t) is obtained as a ground state of ID harmonic oscillator: 

h 2 d 2 Vn moj z 2 z\ 
' n ^T = -^n^ + ^2—^ (22) 

*ll = (^) 1/4 exp [-mu z z 2 /(2h) -iu.t] (23) 
nn 

The transverse part of wavefunction ^±(f±,t) is to be obtained by solving "transverse" GPE: 

3^ h 2 

ot 2m 



4irh 2 a(B) 



m 



/OO rOO 
\*H(z,t)\*dz]/[ \* }l (z,t)\ 2 dz], (24) 
-oo J — oo 



where 
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V± - W ^ |r± ' 2 ~ Re[a(w)] \E(fl)\ 2 (25) 

Because of normalization 

poo poo 

\V\\{z,t)\ A dz = 1/2 and I |*| (z, t)\ 2 dz = 1 (26) 

) J —CO 

The following 2D GPE results from separation of variables for "pancake" trap: 

*fi -sr = -77-A±*± + KL(rl)*L + L^$ ± * ± 2 27 

ot 2m m 

The scattering length a as a function of magnetic field is: 

o(B) - a h3 (1 + (28) 

where As is a width of Feshbach resonance, Br - resonant magnetic field, a ^-background scattering length^. 

Formally the separation of variables is applicable each time when Hamiltonian is a sum of components depending 
on different groups of variables, but this method have additional physical meaning for asymmetric potentials. The 
examples are elongated in z-direction trap and "pancake" 2D trap (fig. [TJ) as in current case. The dynamics of W 
in 2D traps was considered in large amount of papers, including the geometries of periodic potentials, Bessel lattices 
etc. In current case the SFV lattice under consideration has some features, qualitatively described above in discussion 
of classical motion of dipole around phase singularity. 

The continuous transfer of angular momentum from optical vortex to BEC with wavefunction W might induce 
the vortices in initially nonrotating BEC. Because the direction of rotation of classical particle is determined by 
distribution of TC of optical vorices £em of trapping beam the distribution of TC's in quantum superfluid lattice 
(bec will be set identical to those of trapping field. 

The possible solution of equation ([27]) presumes the identical spatial distributions for fields ^±(r±, t) and E(r±). It 
means that ^ is also a sum of a LG functions with alternating topological charges £bec (see eq.pTl)). Unfortunately 
the effective diameter of the core b is much larger than effective size of LG beam bottleneck (4-10 /an). But this 
discrepancy does affect the basic features of solution because of logarythmic dependence of vortex energy on vortex 
dimensions (l7| . 

The correlation K of two complex spatially inhomogeneous fields ^±(f]_,t) and E(r±) is expected to be unit: 



|J^(rl)d'rl|' 
[J\^ ± \ 2 d 2 rl}[J\E(rl)\ 2 d 2 rl] [ ' 



The key point is in adjusting the parameters in ^x(r±,t) and E(f±) in such a way that two last terms in ([24]) 
would cancel each other. This might happen when following condition is imposed upon the coefficients: 

r / , I I o 27rfi 2 a(B) , , 

Re[a(u)} E 2 = ^ 30 

m 

The mutual compensation of these two terms could be achieved via tuning of magnetic field B near Feshbach 
resonance. 

Consider first the case of single vortex trap collocated with a single parabolic subtrap: 

V X = muJ ^ r2 _ Re[a(u)] \E \ 2 r 2 exp [~r 2 /d 2 ] (31) 
The following exact solution for 2D harmonic oscillator is known : 
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which is similar to variational solution for "transverse" wavefunction (|12|) [llj . Note in our case the exact wavefunction 
of the transversal GPE found, rather than variational one p^|) . The stability analysis will be published elsewhere. The 
angular momentum per particle is given by: 

d*± ,2 



<£> = j J *U-ih) -jj^d'rl = h (33) 

where <fi is azymuthal angle (fig. ??). Again SFV carries angular momentum H per particle and kinetic energy of 
the whole vortex line E kin — p irh 2 x ln(b/a)/m [l7[. Consider now the trapping of BEC by of phase-locked Gaussian- 
Laquerre beams placed at the nodes i, j of rectangular grid of period p ( see Eq. |HJ). 

Let us assume for simplicity that optical wavelength A is equal to De-Broiglie wavelength Xdb- 

X = X db = h (34) 
The corresponding BEC temperature for S5 Rb atoms is: 

T = l.7- J K (35) 

After imposing compensation condition (|2"6"]l the residual part of transverse GPE (|2"7|) corresponds to free-space 
evolution: 

ih—± = - — A ± >S> ± + ± tf ± 36 

at 2m 2 

The free space propagation (fig. Q]) equation for EM-field will be of similar form: 

dE i . kr 2 , y 

*=-« A ^ + ^ E (37) 

where f cavity is effective focal length of the laser cavity (fig. [IJ induced by thermal lense collocated with inhomo- 
geneity of optical pumping. 

Because of linearity of equation and superposition principle the solution of ( [37] ) will be the sum of free-space 
modes including zero-order Gaussian functions, Gaussian-Hermit or Laguerre-Gaussian modes located at sites Tj X j v 
[l2| separated by period p. The axes of such optical array are parallel to X,Y. In order to solve numerically eq. 
(|2T|) the split-step FFT method 23] was used. The mesh size in X, Y plane was 512 • 512 points, the "guard bands 
ratio" was chosen equal to 8. So the main part of the field \t was located inside central part of a mesh which had 
64 • 64 size - the "image area". The tolerance of the energy spillover was kept within t\ = 0.0001. The windowing in 
wavenumber space after FFT at each timestep was performed by usage of " Fermi-Dirac" smoothed step function [23j . 
The dissipation inherent to split-step FFT method have led to decrease of total "amount of particles" J J \ ^>\ 2 dxdy 
within "image area" at a speed of 10~ 3 per time step. The special initial conditions of "preselected" SFV array in the 
form (fT7| superimposed upon homogeneous background gave the spatial distribution of " transverse" wavefunction 
well correlated (K sa 0.7) with the OVA array distribution - the array of phase-locked Gaussian-Laquerre wavepackets 
of the first order : 

^ = 2_ [ rn^ L]3/2 E {]f _ rlxJy]) exp ( _ |r1 2/g2) exp { _ rn^-f jx , jy \ » + 

i t B EcArg{r- f jx j y ) + i ir(jx + jy) - i2 uj ± t] (38) 

Thus the numerical solution for ^ j_ of GPE proved to be very close to the linear combination of Gaussian-Laquerre 
functions. Nevertheless this solution takes into account the interference between wavefunctions of " subvortices" , 
because the "doughnut" radius is set to be a bit more than distance between lattice nodes. The arising intereference 
pattern is well correlated with the interference pattern produced by Laguerre-Gaussian OVA |J5(rX)| 2 (fig- E]) with 
the same geometrical parameters and wavelength A. This solution proved to be stable with respect to small harmonic 
perturbations. 



9 



Each SFV carries angular momentum h per particle and rotational kinetic energy Eki n = pTrh 2 \ ln(b/a)/m, b - 
diameter of vortex core, a - interatomic distance, \ - length of vortex which is roughly equal to thickness of "pancake" 
[l7j |. The energy associated with superfluid vortices is of order 10 _ ( 19 "=" 20 ' J or 0.1 ~ 1 eV at density of dilute 85 Rb 
gas p ps 10 16 ^ 18 cm~ 3 . In contrast to superfluid in rotating bucket where angular momenta of vortices are co-directed 
Hi , the BEC vortices trapped by optical vortex array are counter-directed from site to site. Consequently in the net 
sum of angular momenta each vortice of positive topological charge is compensated by the term with negative charge 
and total angular momentum tends to zero. Nevertheless the mutual subtraction of angular momenta(vectorial) of 
adjaicent vortices in the net sum does not mean the mutual subtraction of rotational energies, which are the positive 
scalars. The ground state carry substantial amount of rotational kinetic energy of condensate which contains N 
particles of mass m per unite volume, namely E rot — N vortices pir ft 2 x ln(b/a)/m . 



VI. CONCLUSION 



The optical vortex arrays emitted by solid-state microchip laser were analysed from the point of view of application to 
optical dipole traps. Firstly the numerical modeling of thin slice microchip Fabry-Perot solid-state laser resonator gave 
transverse field distributions (fig. [3]) well correlated with experimentally observed previously. The longitudinal intensity 
distribution consists of periodically spaced array of parallel hollow tubes which slowly diverge while propagating along 
z-axis(fig. [3]). The array of phase- locked Gaussian-Laquerre beams equispaced at the nodes of rectangular lattice 
proved to be a reasonable approximation for experimental and numerical results as well . 

For macroscopic wavefunction \P of BEC trapped in such complex optical field the analytical solutions of Gross- 
Pitaevsky equation were found. These solutions are based upon separation of variables and mutual compensation 
of vortex component of external trapping field via nonlinear term of GPE. The obtained wavefunctions have perfect 
correlation with trapping field, including distribution of topological charges, which form " antiferromagnetic-like" 
lattice. Within framework of this particular model the " antiferromagnetic" lattice of BEC vortices carriers total 
angular momentum close to zero while net rotational kinetic energy of SFV lattice tends to be equal to the sum of 
rotational energies of vortices. Geometrically such BEC - cloud looks like "pancake" perpendicular to z-axis. 

The field of classical velocities, i.e. the field of phase gradients, obtained via Madelung transform, forms "labyrinth" 
structure. It means that trapped atoms move in the "pancake" plane, i.e. x — y plane. The rotation of atom around 
some vortex core is accelerated by EM-torque. The radius of rotation is gradually increased. When particle reaches the 
separatrice it comes to another vortex area. The trajectory of particle in transverse plane (x — y plane) is Mobius-like: 
because the number of optical vortices is finite, in classical picture the particle will return to initial vortex eventually 
after roaming for some time inside trapping " labyrinth" . 

In quantum picture represented via analytic solution of GPE for one trapping vortex and numerical solution for 
optical vortex array the coherent macroscopic wavefunction extends all over OV trapping array with transverse spatial 
dimension of several hundred microns. The complex field of velocities, rotational energy and high degree of correlation 
of SFV wavefunction with OV trapping field promise more resistance to decoherence. 

The qualitative analytic solution supports the basic feature related to trapping of macroscopic particles and to BEC 
trapping: the transfer of OAM from optical field to superfluid. The proposed OVA trap might be interesting from 
the point of view of studies of quantum-classical correspondence. 

The mechanism of imposing the topological charges to BEC vortices by means of manipulating the vorticity of 
trapping optical array could result in demonstration of macroscopic quantum interference phenomena. Evidently there 
are 4 possible topologically equivalent combinations of parameters of solution of OV charges (£em = ±1 ) and 
their relative phases (± it ). Thus there exist 4 wavefunctions \& having an indentical probabability \^f\ 2 distribution 
and different orientation of SFL vortices with respect to physical axes x,y of trapping setup, characterised by £bec 
and their relative phases. The transformation of the one \1/ into another one which have different phase structure is 
equivalent to 90° rotation around z - axis. 
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List of Figure Captions 

Fig.lConceptual view of the near field optical trap. Upper plot: Transverse (in XOY plane) distribution of intensity 
in the near field of solid-state microchip laser ||. Middle: The longitudinal scale extends to 6 Talbot lengths. Z 
- axis is directed along optical axis of microchip laser resonator(bottom) . Additional tightly confined parabolic 
well keeping BEC cloud localized in z - direction is depicted via potential V(z) w (z — z') 2 . Such potential is 
assumed to be superimposed by the other microchip laser beam with cylindrical focusing lens at slightly different 
wavelength of radiation from the range 0.98 4- 2.79/im, in order to avoid interference. Fresnel number has value of 
Nf N vort i ces 2 = 64. 

Fig. 2 Diffractive self-imaging of two-dimensional lattice of 8x8 Gaussian beams with period p — 28/im . The 
longitudinal cross-section at y — plane of intensity distribution I(x, y, z) presented. The lattice is self-reproduced at 
Talbot distance, spatial period division occurs at quarter Talbot distances and central lobe forms outside the Rayleigh 
range. 

Fig. 3 a) Intensity distribution in transverse plane of artificial Laquerre-Gaussian vortex array. Vortices(dark holes) 
are at nodes of rectangular chessboard — like grid. Letters £em denote vortices with alternating topological charges, 
changing sign from one site of lattice to another. b)Longitudinal section of vortex array in the near field in x-z plane 
at y — section. The vortex lines arc parallel, the topological charges £em are flipping from one vortex line to 
another. The wavelength A = lyurra, period p of lattice both in x and y directions is SOfim 

Fig. 4 a)Helicoidal phase surface of Laquerre-Gaussian beam. Atomic dipole having "red" detuning moves along 
phase gradient. Trajectory is located at maximum of intensity. St is the component of Pointing vector tangent to 
helix. The major component of S x is directed parallel to beam propagation, i.e. along Z-axis. b)Hclicoidal phase 
surface of Laquerre-Gaussian beam array. The elementary cell consists of four vortices. The ajaicent vortices have 
arternating topological charges £em and alternating angular momenta. In classcal picture the atomic dipole moves 
along phase gradient of a given vortex next it jumps to another one. 

Fig. 5 a)The argument of macroscopic wavefunction 'J' and corresponding field of velocities obtained via Madelung 
transform . The superfluid vortices form the lattice with alternating topogical charges £bec = ±1- The elementary 
cell consists of four BEC vortices whose locations are identical to the vortices of trapping EM-field. Horizontal pair 
has the same charges £bec which are 7r-shifted with respect to each other, the vertical pair have 7r-shifted —£bec 
charges. The altenating charges £bec make the field of velocities continuous. b) The distribution of the square modulus 
of macroscopic wavefunction ^ obtained via numerical modeling. The location of adjaicent counter-rotating vortices 
£bec shown. 
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FIG. 1: Conceptual view of the near field optical trap. Upper plot: Transverse (in XOY plane) distribution of intensity in the 
near field of solid-state microchip laser [5|. Middle: The longitudinal scale extends to 6 Talbot lengths. Z - axis is directed along 
optical axis of microchip laser resonator(bottom) . Additional tightly confined parabolic well keeping BEC cloud localized in 
z - direction is depicted via potential V(z) w (z — z') 2 . Such potential is assumed to be superimposed by the other microchip 
laser beam with cylindrical focusing lens at slightly different wavelength of radiation from the range 0.98 -r- 2.79/im, in order 



to avoid interference. Fresnel number has value of TV 



64. 
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FIG. 2: Diffractive self-imaging of two-dimensional lattice of 8x8 Gaussian beams with period p = 28^im . The longitudinal 
cross-section at y = plane of intensity distribution I(x,y,z) presented. The lattice is self-reproduced at Talbot distance, 
spatial period division occurs at quarter Talbot distances and central lobe forms outside the Rayleigh range. 
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a) b) 

FIG. 3: a) Intensity distribution in transverse plane of artificial Laquerre-Gaussian vortex array. Vortices(dark holes) are at 
nodes of rectangular chessboard — like grid. Letters £em denote vortices with alternating topological charges, changing sign 
from one site of lattice to another. b)Longitudinal section of vortex array in the near field in x-z plane at y — section. The 
vortex lines are parallel, the topological charges £em are flipping from one vortex line to another. The wavelength A = lfim, 
period p of lattice both in x and y directions is 30/um 




FIG. 4: a)Helicoidal phase surface of Laquerre-Gaussian beam. Atomic dipole having "red" detuning moves along phase 
gradient. Trajectory is located at maximum of intensity. St is the component of Pointing vector tangent to helix. The major 
component of S z is directed parallel to beam propagation, i.e. along Z-axis. b)Helicoidal phase surface of Laquerre-Gaussian 
beam array. The elementary cell consists of four vortices. The ajaicent vortices have arternating topological charges Iem and 
alternating angular momenta. In classcal picture the atomic dipole moves along phase gradient of a given vortex next it jumps 
to another one. 
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FIG. 5: a)The argument of macroscopic wavefunction \& and corresponding field of velocities obtained via Madelung transform 
. The superfluid vortices form the lattice with alternating topogical charges Ibec = ±1- The elementary cell consists of four 
BEC vortices whose locations are identical to the vortices of trapping EM-field. Horizontal pair has the same charges £bec 
which are 7r-shifted with respect to each other, the vertical pair have 7r-shifted —Ibec charges. The altenating charges £bec 
make the field of velocities continuous. b) The distribution of the square modulus of macroscopic wavefunction * obtained via 
numerical modeling. The location of adjaicent counter-rotating vortices Ibec shown. 



